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Abstract 

We construct a model of innovation diffusion that incorporates a spa- 
tial component into a classical imitation-innovation dynamics first in- 
troduced by F. Bass. Relevant for situations where the imitation pro- 
cess explicitly depends on the spatial proximity between agents, the 
resulting nonlinear field dynamics is exactly solvable. As expected 
for nonlinear collective dynamics, the imitation mechanism generates 
spatio-temporal patterns, possessing here the remarkable feature that 
they can be explicitly and analytically discussed. The simplicity of 
the model, its intimate connection with the original Bass' modeling 
framework and the exact transient solutions offer a rather unique the- 
oretical stylized framework to describe how innovation jointly develops 
in space and time. 
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1 Introduction 

Since the middle of the XX th century, a substantial literature em- 
phasizes the importance of quantitative models that enable the fore- 
casting of the diffusion of technological innovations (DTI) not only 
for pure academic interest, but for its practical relevance. While the 
empirical evidence that quantification of DTI is possible was initially 
recognized by E. Mansfield pQ and Z. V. Griliches [2 [3], the first quan- 
titative stylized dynamical modeling framework was proposed by F. 
Bass in his seminal 1969 paper [JJ. Similar in essence to the P.-F. 
Verhulst's epidemiological logistic equation, the Bass' model uses an 
aggregated differential approach enabling the reproduction of the 
relevant dynamical features of the adoption of a new product and/or a 
new technology in a society of consumers. Bass' nonlinear dynamics is 
basically governed by the ratio of two control parameters, namely the 
innovation and imitation rates. Introducing a quadratic nonlinearity 
into the evolution equation, the model mathematically describes con- 
sumer imitative interactions. This nonlinearity leads to an evolution 
characterized by two distinct time scales, i.e. a fast initial exponen- 
tially growing phase followed by a slow asymptotic evolution when 
full equilibrium demand is nearly reached. The seminal works of J. 
A. Schumpeter [5] and subsequently of B. Jovanovic and R. Rob [6] 
support this observation by illustrating the importance of imitation 
waves in DTI and the formation of business cycles. 

In general, interaction-based models have been employed for various 
applications ranging from epidemiology [7], variance of crime rates 
across space when direct interdependencies occur between nearest 
neighbors [8] , herd behavior in financial markets [9[ \TU[ [TT| [T2] to social 
movements and political uprisings (see [13] and |14j for several addi- 
tional references). Interaction-based methods have also been useful 
tools in modeling the diffusion of innovation [TBJ El CE21 EH EH EJJ . In 
their contribution, G. Ellison and D. Fudenberg [22j [23] study agents 
who consider the experiences of their neighbors in deciding which of 
two technologies to adopt, in a world where players use rules of thumb 
that ignore historical data but may incorporate a tendency to use the 
more popular technology. In this learning environment, agents ob- 
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serve both their neighbors' choices, and periodically reevaluate their 
decisions, as opposed to making a once-and-for-all choice. R. Ander- 
gassen et al. |24] investigate the evolutionary process of imitation and 
innovation as a search mechanism in a given neighbourhood of firms. 
In their world, the spreading of information through neighbourhoods 
allows firms to acquire knowledge leading to innovation waves as firms 
attempt to glean information on best practice techniques, which they 
subsequently imitate. For additional models where preference order- 
ings over alternatives in a choice set can depend on the actions chosen 
by other agents, see [251 1261 12T1 128] . 

The original Bass' dynamics is aggregated and hence fails to detail 
the influence of spatial location of agents in the imitative be- 
havior of consumers. Intuitively, spatial considerations strongly affect 
the interactions between agents and the underlying imitation mecha- 
nisms. In fact spatial proximity is argued to be a major driving force 
of innovation diffusion which often exceeds the external marketing ef- 
forts such as advertising |29[ 130] . In 1991, P. Krugman [31] pointed 
out that production is remarkably concentrated in space. This obser- 
vation opened a strong research effort devoted to the understanding of 
the spatial dimension of innovation diffusion. M. P. Feldman [32[ [33] 
echoed P. Krugman's observation in pointing out that geographical 
effects are even more stringent for innovative activity because the ra- 
tionale for the formation of adopter clusters is related to the role of 
word-of- mouth and imitation in the diffusion of innovations. As em- 
phasized in [3l], a clear correlation exists between geographical prox- 
imity and the strength and speed of word-of-mouth spread, sometimes 
labeled as the neighborhood effect. 

In addition to these works, a wealth of empirical studies exemplify 
the importance of geography in the diffusion of knowledge and R & 
D. Spatially- mediated knowledge spillovers of R & D are explicitly 
discussed in [351 Ell EH EE] . It is noteworthy to observe that this pure 
geographic view can be generalized by defining metric distances on 
abstract state spaces in order to describe the evolution of technologi- 
cal advance, R&D investment volume or any other abstract features 
[39] on which agents can compete by adjustment of their individual 
behavior. 

Recent literature suggests that imitation interactions between inter- 
acting agents like bacteria, flies, quadrupeds or fishes can explain the 
formation of compact spatio-temporal patterns, i.e. swarms or pla- 
toons, which spatially evolve as quasi-solid bodies [401 141] . The flock- 
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ing mechanism originates from mimetic type decisions based on agents' 
observations of their neighbors. To the best of our knowledge, spatial 
flocking mechanisms seem to be barely discussed in the interaction- 
based socioeconomic literature. Hence a natural and simple attempt 
to analytically infer the role of spatial parting is to introduce spatial 
effects and it is the aim of our paper to incorporate their influence 
into the original Bass' evolution. 

Adding a spatial dimension transforms the Bass' ordinary differen- 
tial equation into a partial differential equation (PDE). Due to the 
underlying imitation mechanism, the resulting PDE will be intrin- 
sically nonlinear, a perspective that generally offers little hope for 
explicit solutions in the realm of field theories. The present paper 
illustrates how a simple natural spatial extension of the original Bass' 
dynamics leads nevertheless to a fully solvable nonlinear field dynam- 
ics, a truly remarkable result. The resulting equations belong to the 
discrete velocities Boltzmann equations (DVBE) which describe the 
macroscopic properties of a dilute gas. The specific DVBE that can 
be derived from the Bass' dynamics coincide with the Ruijgrok-Wu 
(RW) model introduced and solved by T. W. Ruijgrok and T. T. Wu 
[12]. This intimate connection with statistical physics suggests that 
the Bass' dynamics can be obtained, via a mean-field limit, from a 
microscopic point of view in which a large number of agents interact. 
While the mean-field approach is a basic tool in statistical physics of 
large systems, it has now been explicitly used in recent econometric 
studies as well (see illustrations in QUI QU Q21 O US [28]). Relying 
on the law of large numbers, the mean-field limit allows one to write 
deterministic evolution for probability densities in question. In the 
sequel, we will explicitly construct the microscopic connection that 
exists between the Bass' imitation model with spatial effects and the 
RW model inspired by a similar approach adopted in [43] for a multi- 
agent dynamics in logistics and econophysics contexts. In [33], the 
dynamics exhibit a nonlinear term due to a specific imitation mecha- 
nism giving rise to the famous Burgers' nonlinear PDE to describe the 
emergence of spatio-temporal patterns. As illustrated in [44] . RW dy- 
namics actually generalizes the Burgers' equation, the spatio-temporal 
Bass' model presented here can itself be viewed as a natural general- 
ization of the multi-agent imitation model studied in [43] . 

Besides its direct practical relevance, the simplicity of the original Bass 
model, for which exact analytical solutions are available, has undoubt- 
edly contributed to its popularity in the economics and management 
literatures. Endowing Bass' dynamics with spatially-dependent imita- 



4 



tion mechanisms confers a new dimension to interaction-based socioe- 
conomic modeling, opening the possibility to analytically study the 
generation of spatio-temporal patterns in a highly nonlinear context. 
Our stylized dynamics can be viewed as an exceptional possibility to 
analytically observe the spatio-temporal effects arising for a collection 
of agents subject to imitation interactions. 

2 Spatially-Dependent Imitation Dy- 
namics 

Consider a collection A of N autonomous agents which are in a migra- 
tion process on the one-dimensional real line R. At any time t G R + , 
we assume that the complete population is composed by two types of 
agents A+ and A-, (A = A- U A+), characterized by two associated 
(x, independent migration velocities V+(x,t) and V-(x,t) on R. At 
any time, each agent is free to modify his/her velocity. Agents G A, 
k = 1,2, N, flip their velocities either spontaneously or after an au- 
tonomous decision based on an imitation process (IP). Let us write 
a(x,t), respectively f3(x,t), as the spontaneous transformation rates 
from states A+ ^ a (x,t) <A—, respectively A- >-tp( x ,t) A+- Apart from 
these spontaneous transitions, additional transitions are assumed to 
be triggered by mutual agents' interactions. Specifically, for agent a&, 
the IP mechanism is assumed to depend on the observation of the 
present velocity states, {i.e. V + or V-), adopted by other proximity 
members located in the neighborhood A/& = A4,- LlAfk,+ C A of agent 
afc. For an arbitrary agent dk,±, we define his/her imitation decision 
rule according to his/her interactions with other agents as follows: 

i) Dynamic rule for an agent <2fc ; _. At time t, the agent 
a,k- simultaneously observes the (velocity) states of the agents 
contained in his/her neighborhood A4- The presence of agents 
a,'4- G A/fc, j 7^ k, seen by a^^ triggers an imitation mecha- 
nism which enhances the transition rate towards state V+, i.e. 
P(x,t) i—)- \fl(x, t) + ik(x,t)]. The extra contribution i^{x,t) is 
proportional to (monotonically increasing with) the number 

of agents a Ji+ in the neighborhood A4 of agent cij- (i.e. = 
card (A4,+))- 

ii) Dynamic rule for an agent a^^. At time t, the agent 
Ofc,+ simultaneously observes the (velocity) states of the agents 
contained in his/her neighborhood A/jt and does not modify 
its velocity whatever he/she observes. 
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According to these dynamic rules, the time-dependent position Xk(t) 
of agent ak, k = 1, 2, N, can be written as a set of coupled stochastic 
differential equations (SDEs): 



X k (t) = h{t) 



k = 1,2, ...,N, 



(1) 



where Jfc(t) stands for a two- velocity-states Markov chain (i.e. the 
state space is here O := {V-(x,t),V+(x,t)}) the transition rates of 
which are defined by : 



V+(x,t) ^ a ( X!t ) V-(x,t) V-(x,t) *-t]p( x> t)+i k (x,t)] V+(x,t). (2) 



The noise source in Eq.([TJ) can also be viewed as a non- homogeneous, 
alternating Markov renewal process in which the inverse transition 
rates a(x,t)~ 1 and [j3(x,t) + ik(x,t)]~ l are respectively the mean so- 
journ times in states V + (x,t) and V-(x,t). Observe therefore that 
in Eq.([T]), the coupling between the various agents is realized via the 
extra ik(x,t) transition rate. 

In what follows, we consider a very large population of agents, i.e. 
N — > oo, and we focus on observation neighborhoods Afc, k = 1, 2, N, 
that contain on average = pN agents (0 < p < 1), implying in turn 
that A^. — > oo. Instead of individual trajectories, let us now think 
in terms of statistical properties characterizing the evolution of the 
whole population of agents A. To this aim, we write P(x, t) G [0, 1] 
and Q(x,t) G [0, 1] to denote the density of agents A+ and A- to be 
found at position x G R at time t. 

We now adopt a mean- field approach (MFA), which consists in 
considering that, statistically, the time evolution of an arbitrary agent, 
say ak G A, is representative of the whole population and the influ- 
ence due to his/her neighbors located in A/fe is viewed as an external 
effective interactive mean field. Accordingly, the MFA allows us 
to replace Eq.(pQ) by a single scalar SDE: 



where in Eq.Q, the ik(x, t) transition rates are replaced by an effective 
rate i(x,t) reading as: 



X(t) = I(t) 



(3) 




— — a* + = lim — : a-; + ~ 



pN J\f 




■x+r/2 



P(z,t)dz. 



x-r/2 



(4) 
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where the radius T G IR + characterizes the size of the typical neigh- 
borhood interval M where imitations occur. 



Note that while we proceeded here heuristically, a rigorous mathemat- 
ical justification for the validity of the MFA for this type of dynamics 
relies on the property of propagation of chaos, as explained in [47] for 
SDEs of the type given by Eq.([T]) when the noise is a White Gaussian 
Noise (WGN). SDEs of the type given by Eq.© driven by WGN have 
recently been used for multi-agent dynamics with imitation process in 



Associated with the SDEs given by Eqs.([3|) and (j3|), the MFA en- 
ables us to write a Fokker-Planck equation for the probability densities 
P(x, t) and Q(x, t), 



d 

P(x,t)+V+(x,t)—P(x,t) = +J(P(x,t),Q(x,t))-aP(x,t)+pQ(x,t), 
d 

Q{x,t)+V-(x,t)—Q(x,t) = -J(P(x,t),Q{x,t))+aP(x,t)-PQ(x,t), 
ox 

(5) 

where the imitation rate term in Eqs.© reads as: 

rx+F/2 

P(z,t)dz . (6) 



J(P(x,t),Q{x,t)) = Q(x,t) 



x-T/2 



The dynamics in Eqs.([5]) is a coupled set of nonlocal and nonlinear 
field equations, which barely offers hope for any analytical discus- 
sion. However, for small p implying small radius V {i.e. infinitesimal 
interaction neighborhoods), we may Taylor expand Eqs.(0), up to first 
order in T, to obtain: 



P(x,t)+V + (x,t)—P(x,t) = +TP(x,t)Q(x,t)-aP(x,t)+f3Q(x,t), 
ox 

d 

Q(x, t) + V-(x, t)—Q(x, t) = -TP(x, t)Q(x, t) + aP(x, t) - f3Q(x, t). 
ox 

(7) 

We directly observe that Eqs.([7]) can be viewed as a generalized Bass' 
dynamics which confers relevance to the spatial dimension on which 
agents evolve. Indeed, similarly to the original Bass' model, we in- 
clude an imitation process, represented in Eqs.([7|) by the nonlinear 
contribution TP(x, t)Q(x, t). 
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Writing P(x,t) + Q(x,t) = the summation of both equations 

in d?]) yields a continuity equation: 

^(x,t) + [V + (x,t) + V.(x,t)} ^(x,t) = 0, (8) 
which in turn implies the normalization constraint: 

/ P(x,t)dx+ / Q(x,t)dx = l, Vt€M + . (9) 

Although the dynamics given by Eqs.(J7J) has been derived for non- 
homogeneous and non-stationary parameters a(x, t), /3(x, t) and V±(x, t), 
in the sequel we will restrict our attention to situations where these 
parameters can be assimilated to constants (i.e. a, (3 and V±). 

2.1 Spatial Homogeneous Regimes - Bass' Model 

When V- = V+ = 0, the spatial character disappears from the dy- 
namics given by Eqs.©, which implyies that P(x,t) = P(t) and 
(Q(x,t) = Q(t). More precisely, Eqs.([7]) becomes: 

P(t) = +TP(t)Q(t) - aP(t) + /3Q(t), 

Q(t) = -TP(t)Q(t) + aP(t) - (3Q(t), (10) 

with the notation P(t) := ^P(t). The constraint given by Eq.([9]) now 
simply becomes P(t) + Q(t) = 1 and enables us to rewrite Eq. ffTO]) in 
the following form: 



p(t) = -r\p 2 (t) 



1 (a + P) 

r 

At this stage, it is worth observing that for a = and T = 1, Eq. tjlip 
reduces to: 

P(t) = [1 - P(t)] I? + P(t)] , (12) 

which is precisely the original Bass' dynamics [1] with /3 being the 
ratio between the imitation and innovation rates. 



For completeness of the exposition, let us integrate Eq. ljlip . with the 
initial condition P(t = 0) = Pq, to get: 



pm = (A + b)(P -A + b)e- 2Art + (A-b)(P + A + b) 
[ > (P + A + b) - (P - A + b) e" 2Ar * { 6) 

with the definitions: 



S 



1 _ . . 1 /,„ „, 2 



b = -—(T-a-(3) and A = —y (V — a — 0) + 4/3T. 
In the asymptotic time limit t — > oo, Eq. (fT3]) converges to: 

lim P(t) = P station = A - b. (14) 



t— >oo 



For the original Bass' model obtained when a = and T = 1, we have 
b = -|(1 - /3) and A = ±(1 + /?)• Accordingly, with P = 0, Eq.([T3]) 
reduces to the original Bass' solution: 

Pit) ~ /? [ 1 - £ " (1+/3)t ] Q(t) ~ ( 1 + ^~ (1+/3) * (15) 

Moreover P s tation = A — 6 = 1, which expresses the fact that all agents 
ultimately adopt the new technology, as it is obviously expected in 
the original Bass' modeling framework, when a = 0. 



3 Bass' Dynamics with Spatio- Temporal 
Effects 

Coming back to the general dynamics given by Eqs.([7|) and introduc- 
ing dimensionless coordinates via the Galileo transformation (x, t) \— > 
(y, s) defined by: 




(16) 



we straightforwardly have: 



d x {-) i— >■ 



2r 



(V+-V-) 



dy(-) and d t {-) h-> 



(V- - V+) 



d y (-)+Td s (-), 
(17) 

which transforms Eqs.([7]) into the coupled set of nonlinear partial dif- 
ferential equations (PDEs): 

a 

8 S P + d y P = +PQ ~fP + ^Q, 



d s Q - d y Q = -PQ + ^P-^Q 



(18) 
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The set of non-linear PDEs given by Eqs. (|18p can be interpreted as 
being a discrete two velocity model of Boltzmann equations, first stud- 
ied in |42j . The dynamics given by Eqs. (|18p is remarkable, as using a 
generalized Hopf-Cole logarithmic transformation, we get: 

P(y, = + d s log H(y, s) - d y log H(y, s) (19) 

Ol 

Q(y, s) = --d s log H(y, s) - d y log H(y, s) (20) 
which actually reduces the set Eqs. (fT8j) into the Telegraphist equation: 



(21) 



a/3 

d ss H(y, s) - d yy H(y, s) - -^2 H (y^ s ) = °- 



Accordingly, the dynamics given by Eqs.0 with constant parameters 
has the 1 property that it can be exactly solved for any initial 
conditions Po(y) and Qo(y). According to [42j . the general solution 
reads as: 



H(y, s) = ± [A(y + s) + A(y - s)} + ^(y, s) + ^ B 2 (y, s) 



(22) 



where rj = ^IsL anc [ where we have the following definitions: 

ry+s 



Bi(y,s) 



h (vV* 2 ~(y- y') 2 ) B(y')dy', 



and 



B 2 {y,s) 



y+s 



1 



h(vVs 2 -(y-y') 2 )A(y')dy', 



'y-s \V s2 ~(y- y') 2 / 

with Io( - ) and Ii(-) being the modified Bessel's functions and: 
1 



B(y) 



Po(y) - Qo(y) + £ + £ 



My), 



(23) 



A(y) = exp 



1 /■» 



Po(y') + Q (y')-^ + ^ 



dy' 



(24) 
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3.1 Behavior of the Solutions 



Though fully explicit and exact, the solution given by Eqs.(|19p and 
(|20p deserves discussion and interpretation for specific situations and 
this is precisely the objective of this section. In what follows, we will 
systematically choose V = 1 and unit velocities V± = ±1. 

Let us here focus on the Bass' dynamics that results when a = 0. 
This directly implies that 77 = y/a/3 = and hence Eqs. (|19p and (|20|) 
coincide with the ordinary wave equation. By definition of Bessel's 
functions [H], we have that Io(0) = 1 and li(0) = 0, thus leading to 
the usual wave solution in the form: 

H(y, s) = i [A(y + s) + A(y ~s)} + ^ J" ' B{y')dy> '. (25) 

Due to the fact that a = 0, we expect that the agents' population A+ 
with velocity = +1 increases by opposition to the population A- 
of agents with velocity V- = — 1 which is doomed to extinction. Let 
us now explore the transient nature of the solution and this for three 
types of initial conditions. 

a) Initial conditions: Po(y) = S(y) and Qo(y) = 0. 

In this case, the time-dependent solution reads as: 

P(y,s) = 5(y-s) and Q(y,s) = 0. (26) 

as it can be checked directly from Eqs. ([T9j) and ([20j) . Eqs. ([26"P 
are clearly consistent with the fact that a = and hence that no 
transitions from V+ = +1 to V~ = —1 velocities occur. Hence, 
starting with all agents with velocity V+ = +1, they stay with 
their original velocity and the density P(y, s) is a uniformly trav- 
eling Dirac mass with velocity = +1 towards the positive 
IR-axis. 

b) Initial conditions: Po(y) = and Qo(y) = S(y). 

Either by direct substitution into Eqs. (|18p or alternatively by 
using Eqs. (|19p and (|20p . one can verify that the time-dependent 
solution in this case reads as: 

P(y,s) = ^y-^e(\y\-s) and Q(y, s) = e^ s S(y+s) (27) 

where @(\y\ — s) is a Heaviside cutoff function which identically 
vanishes for negative arguments. The dynamics P(y, s) of the 
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Figure 1: Spatial dispersion P(y,s) of the agents with velocity V + = +1 
for different times s = [0.01; 20; 40]. The initial conditions are given by 
Pq{v) — and Qo(y) = 8(y), a = and j3 = 0.1. In this case, the agents are 
immediately segregated and there are no imitation processes. 



spatial dispersion of the agents with velocity V+ = +1 is illus- 
trated in Fig. [TJ 



The behavior of the solution given by Eqs. (|27p can be easily un- 
derstood. Indeed, the Dirac mass for Q(y,s) expresses the fact 
that agents with velocity V_ = — 1 are gradually depopulated 
at the rate s/3 for the benefit of agents traveling with velocity 
V+ = +1 and hence migrating to P(y, s). The Heaviside function 
&(\y\ — s) expresses the fact that no agent can possibly be found 
at a distance larger than \y\ at time t (remember that the veloci- 
ties here are V± = ±1). It is worth observing that in the original 
Bass' model, the adoption rate given in Eq. (|15p is equal to (1+/3). 
This overcomes the f3 adoption rate that we will obtain below 
in Eqs. (|28p . due to the fact that in the present configuration, 
imitation process does not enter into play in our spatial model, 
which drastically moderates the overall adoption rate. Indeed, 
in this limiting regime the agents' populations ^4+ and are 
immediately and definitively segregated, which hence never al- 
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lows imitation processes to take effect (the unit rate discrepancy 
between the (1 + /3) adoption rate occurring in Eq. (fT5|) and the f3 
rate in Eq. (|27p is obtained from the T = 1 choice). The resulting 
temporal evolution of the agents' overall adoption rate obtained 
in the present case (P(t) = L P(y, t)dy), compared with the one 
of the original Bass' model, will be illustrated later in Fig. |U 



Finally, let us observe that, from Eqs. (|2T|) . one immediately 
obtains: 

/ P(y, s)dy =\l- e"H and / Q(y, s)dy = e^ s , (28) 

thus showing that Eq.([9]) is fulfilled. In addition, for asymptotic 
times s — > oo, Eqs.([28]) indicate that all agents ultimately adopt 
the V+ = +1 velocity as it is expected for the a = regime. 

c) Initial conditions: Po(y) = Qa{y) = ^[tanh(y + 7) — 
tanh(y - 7)] . 

For these initial conditions and for the particular choice (3 = 2 
and 7 = 4, the resulting time-dependent solution of our spatial 
Bass' model is given by: 

^ coshir^+T) (i " tanHy ~ s ~ i)+ tanHy - s + 7))] } 

(29) 

and 

Q ^ s) = l{Hh) [ ~ A(y+s) + 

e „^ cosh(j/ + S -7) (1 _ tanh(y + g _ l) + tanh(y + s + 7))] j 
cosn(y + s + 7) J 

(30) 

where 

H(y, s) = i [A(y + s) + A(y - s) + Z?i(y, a)] , 
A(y) = e y exp / ^ [tanh (z — 7) — tanh (z + 7)] <iz | : 



, cosh(y - 7) 
cosh(y + 7) 



e 



-7 



B±(y, s) = — ^— [atan (sinh(y + s + 7)) — atan (smh(y — s + 7))] 

13 



-e 37 [e-y- 3 ^ + atan(e J/+s+7 ) - e"^" 7 - atan(e s/ - s+7 )] . 

The spatio-temporal dynamics P(y, s) of the agents with veloc- 
ity V+ = +1 is illustrated in Fig. [2j The joint evolutions of 
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Figure 2: Spatial dispersion P(y, s) of the agents with velocity V + = +1 
for different times s = [0; 0.2; 3; 10]. The initial conditions are given by 
Po(y) = Qo(y) = ^ [tanh(y + 7) - tanh(y - 7)], a = 0, /3 = 2 and 7 = \. 

P(y,s) and Q(y,s) (representing agents with velocity V+ = +1 
and V- = — 1 respectively) are drawn in Fig. [3l 

In the present configuration, the two types of agents have an 
identical initial spatial distribution (i.e. Po(y) = Qo(y), Vy). 
Hence, half of the agents have initially the velocity V+ = +1, 
the other half having the velocity V- = — 1 (i.e. f K Po(y)dy = 
f R Qo(y)dy = |). The spatio-temporal behavior of the solu- 
tion given by Eqs. (|29p and (|5D|) can be split into two different 
time phases. For short times of the dynamics, the overlap be- 
tween P(y,s) and Q(y,s) is non-null (i.e. the two populations 
of agents ^4+ and A- are not spatially segregated), thus lead- 
ing to strong imitation processes between the agents. Agents 
are changing their velocity from V- = — 1 to V+ = +1 due to 
both imitation and spontaneous transitions. Accordingly, Q(y, s) 
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Figure 3: Spatial dispersions P(y, s) and Q(y, s) of the agents with velocity 
V + = +1 and V- = —1 respectively, for times s = [0;0.3; 1.2]. The initial 
conditions are given by Po(y) = Qo(y) = §- [tanh(y + 7) — tanh(y — 7)], 
a = 0, j3 = 2 and 7 = |. 

is gradually depopulated for the benefit of P(y,s). The imita- 
tion processes are decreasing as the overlap between s) and 
Q(y, s) gets smaller, but the rate f3 of spontaneous transitions 
from V- = — 1 to V+ = +1 remains itself constant with time. In 
the time asymptotic regime, almost all the agents have adopted 
velocity V+ = +1, leading P(y,s) to behave as a probability 
density uniformly traveling with velocity V+ = +1 towards the 
positive M-axis. 

In Fig. HI the temporal evolution of the overall adoption rate 
P(t) = L P(y, t)dy of the agents is illustrated and compared to 
the one observed for the Bass' original model. For the present 
configuration, the overall adoption rate stands between f3, the 
rate observed for our spatial Bass' model in absence of imitation 
processes (Section 3.1.b), and ,9 + 1, the rate obtained for the 
original aggregated Bass' model. Hence, in general, the overall 
adoption rate of our spatial Bass' model will be equal to /3 + e(t), 
e(t) e]0, 1], with e(t) depending on the initial distributions of the 
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Figure 4: Temporal evolution of the overall adoption rate P(t) i) for the 
original Bass' model, ii) for our spatial Bass' model when there are imitation 
processes (Section 3.1.c) and Hi) for our spatial Bass' model when there are 
no imitation processes because the two populations of agents A+ and A- are 
immediately segregated (Section 3.1.b). 

agents and on the model parameters. Remember that the im- 
itation rate (and hence the overall adoption rate) is controlled 
by the number of neighbors that each agent effectively observes 
during the imitation process. In the aggregated (original) Bass' 
model, this number is maximum as each agent systematically 
observes the global population of agents. In the limiting regime 
of Section 3.1.b, the number of observed agents is equal to as 
the two populations of agents A+ and A- are initially and hence 
permanently segregated, thus allowing no imitation processes. 

4 Conclusion and Perspectives 

Often, agents' remoteness may naturally reduce the efficiency of im- 
itation processes, a feature that is totally absent in the original ag- 
gregated Bass' approach. Incorporating a spatial dimension into the 
Bass' dynamics is however not a minor extension. It transforms in- 
deed a nonlinear single dimension dynamics into a nonlinear infinite 
dimensional field dynamics for which, in general, no solution methods 
are available. It is hence remarkable that our spatial generalization 
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of the original Bass' dynamics leads to a class of models for which 
the evolution of the resulting spatio-temporal patterns can be exactly 
calculated. Indeed, the quadratic nonlinearity, due here to the under- 
lying imitation mechanism, coincides with the collision term found in 
a solvable Boltzmann equation, the RW dynamics, used for gas models 
in mathematical physics. Our paper points out that the RW dynamics, 
solvable via a linearizing logarithmic transformation, offers a unique, 
synthetic and exact modeling framework to study nonlinear features 
generated by spatio-temporal imitation processes in economics sys- 
tems. 

In this paper, we mainly focused on Bass' dynamics which does not 
allow back transitions (i.e. a = 0, implying no transitions from V+ to 
V-). Allowing a > 0, regimes involving shock waves with propagating 
velocities w emerge and the velocity range is < w < 1, [12] . When 
a dominant spontaneous tendency to stay in the V- state exists (i.e. 
a > f3), imitations enhance the spontaneous /3-flow from V- i— > V+ 
and may lead to a time-independent (i.e. w = 0), shock type inhomo- 
geneous solution. Depending on the initial conditions, this marginal 
stationary w = solution separates two regimes of shocks propagating 
with either positive or negative velocities, a rich dynamical behavior 
that deserves further investigations in economy. 
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